3-D hybrid kinetic modeling of the interaction between 
the solar wind and lunar-like exospheric pickup ions in 
case of oblique/quasi-parallel/parallel upstream 
magnetic field 

A.S. Lipatov * 1 ’* 

GPHI UMBC/NASA GSFC, Code 673, Greenbelt, MD 20771, USA 
and Faculty of Problems of Physics and Power Engineering, 

Moscow Institute of Physics and Technology, Russia 

W.M. Farrell, J.F. Cooper, E.C. Sittler Jr., R.E. Hartle 2 
NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA 


Abstract 

The interactions between the solar wind and Moon-sized objects are determined 
by a set of the solar wind parameters and plasma environment of the space 
objects. The orientation of upstream magnetic field is one of the key factors 
which determines the formation and structure of bow shock wave/Mach cone 
or Alfven wing near the obstacle. The study of effects of the direction of the 
upstream magnetic field on lunar-like plasma environment is the main subject of 
our investigation in this paper. Photoionization, electron-impact ionization and 
charge exchange are included in our hybrid model. The computational model 
includes the self-consistent dynamics of the light (H + , He + ) and heavy (Na + ) 
pickup ions. The lunar interior is considered as a weakly conducting body. Our 
previous 2013 lunar work, as reported in this journal, found formation of a triple 
structure of the Mach cone near the Moon in the case of perpendicular upstream 
magnetic field. Further advances in modeling now reveal the presence of strong 
wave activity in the upstream solar wind and plasma wake in the cases of quasi- 
parallel and parallel upstream magnetic fields. However, little wave activity is 
found for the opposite case with a perpendicular upstream magnetic field. The 
modeling does not show a formation of the Mach cone in the case of 6b,u ~ 0°. 
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1. Introduction 

The solar wind plasma, lunar surface regolith, and neutral exosphere are all 
intimately interconnected in complicated ways. For example, solar energy can 
release neutrals from the surface, but these neutrals can become ionized and pro- 
vide photoions that feed back and affect the incoming solar wind. In this paper, 
we continue to examine this newly appreciated environmental system feedback 
pathway. Earlier observations, as reviewed in Lipatov et al. (2013b), showed the 
existence of neutrals such as He, Na, K, and O. Detected ion species included 
H + , He ++ , He + , C + , 0 + , Na + , K + , and A r + . The experiments on board of the 
Lunar Atmosphere and Dust Environment Explorer (LADEE) also provide new 
information about He, A r, Ne, potassium, sodium, OH, H20, Si, Al, Mg, Ca, 
Ti, and Fe neutral components (see e.g. Benna et al. (2014), Colaprete et al. 
(2014)). 

There is a long history of plasma simulations of the lunar environment, as 
also reviewed in Lipatov et al. (2013b), the next step in our work being to 
include photoion pickup and reflected ions and consider how these substantial 
populations affect the solar wind flow. While magnetohydro dynamic (MHD) 
models have been useful for a study of the global configuration of the plasma 
structures near the Moon, important kinetic effects have not been included in the 
MHD formalism: namely, anisotropy of the ion velocity distribution resulting in 
excitation of the low-frequency electromagnetic waves, formation of the electron 
and ion beams and excitation of the high-frequency waves, etc. These effects 
may be recovered by using hybrid or full kinetic modeling. These models take 
into account finite ion gyroradius effects and provide an understanding not only 
on the solar wind flow past the object, but also of the solar wind- pickup ion 
interaction. 

We apply various orientations of the incoming magnetic field ( 6b,u = 60°; 
45°; 0°) to a model of the neutral exosphere with H, He, Na species. We use 
a code which combines a Boltzmann’s “particle-in-cell” approach (Lipatov et 
al., 1998), together with a hybrid plasma model (Lipatov et al., 2002b) in three 
spatial dimensions (see, e.g. Lipatov and Combi (2006)). Charge exchange 
between upstream and pickup ions and the exospheric neutrals is also included. 


2. Hybrid model description 

A quasi- neutral hybrid model is used for a kinetic description of the upstream 
and implanted pickup ions, and a fluid description of electrons, to study the 
interaction between the solar wind and the ionized and neutral components 
of the lunar environment (Lipatov, 2002; Lipatov et al., 2013b). This model 
describes well the wave-particle interactions on the ion spatial and time scales 


2 


(pci = Uq/£1\ and cj < Hi), where p cl is the gyroradius for ions. Uo and fii 
denote the upstream velocity and the ion gyrofrequency. Our model takes into 
account ionization and charge exchange. The finite conductivity of the Moon’s 
interior is also included. 

Figure 1 presents an illustration of the interaction between the solar wind 
and the Moon. In our modeling coordinate system, the X axis is directed away 
from the Sun, Y axis is oriented in the direction of Earth’s orbital plane and Z 
axis completes the right-handed system. Note that the position of the center of 
the Moon is x = 0, y = 0, z = 0. 

Mass and charge state of the ions under consideration are M s and Z s — 1, 
where the ion population subscript s is s = 1 for H + upstream ions and s = 2,3,4 
for H + , He + , and Na + pickup ions. Charge exchange between ions and neutral 
particles is computed with Eqs. (12-15) from Lipatov and Combi (2006). 

In the following we review for convenience the electron and electric field 
formulations from Lipatov et al. (2013b). The electric field is computed from a 
standard generalized Ohm’s law (Braginskii, 1965): 

E = _VPe + JeXB^ (1) 

en e en e c 

The right-hand side of this equation consists of terms with the electron pressure 
p e = nm e (v'^) / 3 = n e T e and the electron current J e , where v' e denotes the 
electron thermal velocity. 

The electron density is computed from the quasi-neutrality condition 


n e — n H + )Up + n H + ?PI + n He +,pi + ^Na+,Pi- 


( 2 ) 


We also apply an adiabatic process for the partial pressures of all electron 
populations: 

/3e,PI,s 5/3 


Pe oc (nf/p + 




), 


(3) 


where /3 e and /? e ,pi, s denote electron betas for the upstream and pickup popula- 
tions, respectively. We also assume here that n e , u P = ni jUp , and n e ,pi, s = ?h,Pi,s- 
The ion kinetic approach accurately includes effects of mass loading, ion pres- 
sure anisotropy, exospheric penetration by ions, and asymmetry of lunar plasma 
flows. 

The electron stimulated desorption (Mclain et ah, 2011) and photoionization 
are assumed to be the main mechanisms for ion production from the neutral 
exosphere of the Moon. 

The equation for production of new ions from a three-species neutral exo- 
sphere near the Moon is similar to our four-species exosphere model (Lipatov 
et ah, 2013b) and corresponds to 


-fexo,k — ^ / z,/c^ / exo,k^ 7(^5 0) GXp [(r ^exobase)/^exo,k] • (4) 


Here, n exo ,/c is the maximum value of the neutral density at the exobase, where 
v exobase ~ i?M (lunar radius). Index k denotes exospheric species. The ex- 
ospheric densities, spatial scales, initial thermal velocities of pickup ions and 
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effective ionization rates are chosen in the same manner as in Lipatov et al. 
(2013b). The maximum values of the neutral density at the exobase are n exo ,H = 
1.0 x 10 2 cm -3 , n exo ,He = 2.0 x 10 3 cm -3 , and n exo ,Na = 1.0 x 10 2 cm -3 . Spatial 
scales are: h e xo,Na = 87 km, h exo ,H = 1990 km, h e xo,He = 498 km (Hartle and 
Killen, 2006) and v t h, P i = 2.0km/s. 

Effective values of the ionization per atom or molecule are: v x H = 7.5 x 
10 -7 s -1 , z/i 5 H e = 7.8 x 10 -9 s _1 (Wang et ah, 2011b), and zq ? Na = 1.8 x 10 -5 s _1 
(Potter and Morgan, 1988; Gruntman, 1996; Wang et al., 2011b). Total pickup 
ion production rates are: G H + ~ 5.0 x 10 21 s _1 , G He + ~ 1.0 x 10 21 s _1 , and 
Gn & + ~ 2.0 x 10 22 s -1 . These values may be converted to the following mass 
production rates: mn+Gn+ ~ 8.35 x 10 3 g/s, m He +G He + « 6.68 x 10 3 g/s, 
and ?«Na+ G]sr a + ~ 7.682 x 10 ~ 1 g/s. The pickup ions of variable weight are 
generated with algorithm described in Lipatov (2012). 

At the flank boundaries a damping boundary condition for the electromag- 
netic field has been applied. At the downstream boundary magnetic field pertur- 
bations and particles may escape the computational domain. A portion of the 
particles may re-enter into the computational domain from the outflow plasma. 
The upstream values of the magnetic field and electric fields are the same as in 
the solar wind. 

We assume that ion bulk velocity equals to zero inside the Moon. The finite 
conductivity of the lunar core is taken into account with the effective magnetic 
Reynolds number R m ,core = 47rGoTcr e ff, C ore/c 2 . 

When pickup ions pass through the surfaces y = ±(0.5DY — 5 x Ay), or 
z = ±(0.bDZ— 5x Az), or x = 16L — 5 x Ax we remove them from computational 
domain. The particles are absorbed when they intersect the lunar surface. We 
do not use any boundary condition for the electromagnetic field at the lunar 
surface. Inside the Moon, we use our equations for the electromagnetic field 
with internal conductivity and the ion bulk velocity. The variations of the value 
of the conductivity and ion bulk velocity across the lunar surface create a jump 
in the electric field. 

A computational domain is chosen with following dimensions DX = 21 L, 
DY = 30 L, and DZ = 20 L. Here the spatial scale L equals the lunar 
radius Rm • The computational domain contains a grid of 211 x 301 x 201 
points. 10 9 macro-protons and 10 8 macro-pickup ions were used in the mod- 
eling. The numerical time step At is chosen to satisfy the condition At < 
min(Ax, Ay, Az) / (16v max ) , where r TOax denotes the maximum value of the 
macro-particle velocity. 

A homogeneous solar wind plasma flow was used to start a modeling. An 
external electromagnetic field has upstream value: B = Bo and E = — Uo x 
Bo/c. In the lunar interior the electric field is E =0 and the magnetic field is 
B = Bo- The ion bulk velocities also equal to zero. 

For chosen values of the upstream magnetic field and the local bulk velocity 
the values of the gyroradius of pickup ions (He + , Na + ) and upstream ions (H+) 
are about 800km, 4500km and 200km (respectively). According the analysis 
from Winske et al. (1985) we have good resolution for the low-frequency waves 
and insufficient resolution for the high-frequency waves. 
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3. Results 


In our simulation we have adopted the following input parameters. The solar 
wind velocity is Uo = 305km/s, U = (1.0, 0.0, 0.0) Uo] the density is no = n H + = 

3.0 cm -3 and the value of the magnetic field is Bo = 5.2 nT. In oblique case 

B = (—2.6, 0.0, 4.5) nT; in quasi-parallel case B = (—3.69, 0.0, 3.69) nT and 
in parallel case B = (—5.2, 0.0, 0.0) nT. The Alfven and sonic Mach numbers 
are: = 5.2 and Ms = 3.7. Plasma betas are: /3 H + = 0.2, /3 e = 0.5 and 

/^pickup, e = 0.05. Inside the undersurface shell the effective Reynolds number is 
Rm, shell = 0-05 whereas inside the lunar conducting core is R m ,core = 0-2. Here, 
the lunar conducting core has a radius of Rm/2. 

Our modeling has been performed up to ten average transition times for par- 
ticles from the upstream boundary to the downstream boundary, t = 10T trans i t . 

3.1 Solar wind-Moon interaction in case with 0 u,b = 00° 

Consider first a case for oblique upstream magnetic field. Figure 2 (top) 
shows the cuts of Na + pickup ion density. The heavy Na + pickup ions create an 
asymmetrical plasma wake due to rotation around the external magnetic field. 

Figure 2 (bottom) shows the density of H + pickup ions. The modeling shows 
a formation of asymmetrical halo with light pickup ions near the Moon. The 
plasma wake has an orientation oblique to the incoming magnetic field. The dy- 
namics of light pickup ions are not similar to that observed in the perpendicular 
interaction Lipatov et al. (2013b). 

Figure 3 (top) shows the solar wind H + ion density. The background plasma 
creates an asymmetrical plasma wake due to mass loading by the heavy pickup 
ions creating increased ion density on the outer boundary (y < 0, Fig. 3, top). 
The central region of the plasma wake is shifted slightly in both y and z com- 
pared to the nominal solar wind flow. 

At the exobase the solar wind density is nn+ )max ~ 1.0 no in the case with 
an oblique magnetic field. The density of various species of pickup ions is: 
n H + ~ 2.0 x 10 -4 no, n He + ~ (1.5 — 3.0) x 10 -5 no, n Na + ~ 3.0 x 10 -4 no- 
Normalization density is no = 3 cm -3 . A jump in the proton density across the 
Mach cone is about 0.3 no inside the cut located at x/L = 3.0. At the larger 
distance (x/L > 10) the jump in proton density becomes weaker. 

The magnetic field (B z ) profile has also an asymmetrical configuration with 
reduced value in upper part (y > 0) and increased value in lower part (y < 0), 
(Fig. 3, bottom). In the day side of the exosphere one can see a jump in 
the value of the magnetic field (magnetic barrier), A B z ^ 0.2 Ho. A jump in 
the magnetic field across the Mach cone is about 0.3 Ho- There are no visible 
fluctuations in the magnetic field in upstream and plasma wake. 

The feature of our modeling of the solar wind-Moon interaction in case of 
oblique upstream magnetic field is the formation of a strong asymmetrical Mach 
cone, asymmetrical heavy and light ion plasma wakes, and asymmetrical induc- 
tive magnetotail. The mass loading of exo-ions should have a real and observable 
effect on the wake plasma and magnetic structure. Onboard measurements of 
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ions and magnetic field in the lunar plasma wake will allow to make a conclu- 
sion about the composition of pickup ions and their sources. This modeling 
prediction makes this application unique. 

3.2 Solar wind-Moon interaction in case with Qu,b = 45° 

Consider now the case of a quasi-parallel upstream magnetic field. Figure 
4 demonstrates 2-D cuts of the (H + , Na + ) pickup ion density. The x — y cuts 
(left column) are located at z = 0, y — z cuts (middle column) are located at 
x/L = 3.0, and x — z cuts (right column) are located at y = 0. Let us consider 
first the distribution of light pickup ions. Figure 4 (bottom) shows the cuts of 
the distribution of the H + pickup ions density. Pickup ions create an asym- 
metrical halo around the Moon. The plasma wake has an asymmetrical form 
and an oblique orientation along the incoming magnetic field. This asymmetry 
is created by the polarization electric field in the exosphere. A more diffusive 
structure compared to the oblique case may be explained by formation of the 
shell velocity distribution and wave activity in this case. 

The Na + pickup ion density is shown in Figure 4 (top). The x — y cuts 
(left column) are located at z = 0, y — z cuts (middle column) are located at 
x/L = 3.0, and x — z cuts (right column) are located at y = 0. The Na + pickup 
ion creates an asymmetrical wake due to gyroradius. Figure 5 (top) shows the 
solar wind H + ion density. The background plasma wake has an asymmetrical 
plasma wake due to mass loading by the heavy pickup ions with increased ion 
density on the outer boundary (y < 0, Fig. 5, top). The central region of the 
plasma wake is mildly oblique to the solar wind bulk velocity. 

Figures 4 and 5 gives the following variation in the densities of background 
and pickup ions at the exobase in case with quasi-parallel magnetic field: n H + ~ 
1. no, n H + « 1.2 x 10 -4 no, n He + « (1.5 — 3.0) x 10 -5 no, n Na + « 2. x 10 -4 no- 
no = 3.0 cm -3 is upstream solar wind density. A jump in the proton density 
across the Mach cone is about 0.15 no- 

The magnetic field (B z ) has reduced value in upper part (y > 0) and in- 
creased value in lower part (y < 0), (Fig. 5, bottom). A jump in the magnetic 
field across the Mach cone is about 0.15 Bo, and it is smaller than in the case 
of @b,u = 60°. Upstream, 1-D cuts show the following perturbations in the 
magnetic field with amplitude about of SB y ~ 0.01F?o and 5B Z « 0.01i?o, and 
wave-length about of A ~ 1400 km. Downstream, 1-D cuts show much stronger 
perturbations in the magnetic field with amplitude about of SB y « 0.01F>o and 
5B Z ~ 0.01F>o, and wave-length about of A ^ 1200 km. 

Once again, we conclude that for a quasi-parallel B-field configuration, the 
mass loading of the exo-ions creates a modified solar wind flow, changes the 
orientation of the plasma wake and alters the structure of the tailing B-field 
(compared to the ’unloaded’ situation). Pickup ions have a clear and obvious 
effect on the near-Moon plasma interaction. 

3.3 Solar wind-Moon interaction in case with Qu,b ~ 0° 

Consider now the case of a parallel upstream magnetic field. Figure 6 demon- 
strates 2-D cuts of the (H + , Na + ) pickup ion density. The x — y cuts (left col- 
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limn) are located at z = 0, y — z cuts (middle column) are located at x/L = 3.0, 
and x — z cuts (right column) are located at y = 0. Let us consider at first 
the distribution of light pickup ions. Figure 6 (bottom) shows the cuts of the 
distribution of the H + pickup ions. Pickup ions create a halo around the Moon 
which is transformed into the relatively thin tail along the lunar plasma wake. 

The Na + pickup ion density is shown in Figure 6 (top). The x — y cuts 
(left column) are located at z = 0, y — z cuts (middle column) are located at 
x/L = 3.0, and x — z cuts (right column) are located at y = 0. The Na + heavy 
pickup ions form a halo near the Moon and a complex tail structure rotating 
around the magnetic field. The trajectory of the heavy Na + pickup ions in 
the outer tail is similar to a test particle motion with an average value of the 
gyroradius about of 3 x 10 4 km. 

The modeling does not demonstrate a formation of the Mach cone in the 
solar wind H + ion density profile (Fig. 7, top). The central region of the plasma 
wake is slightly oblique to the solar wind bulk velocity. 

Figure 7 (bottom) shows a distribution of the B z magnetic field. This figure 
demonstrates an asymmetrical tail and strong perturbations in the magnetic 
field but without a formation of Mach cone. Upstream, 1-D cuts show the 
following perturbations in the magnetic field with amplitude about of SB y « 
0.02F>o and SB Z ~ 0.02F>o, and wave-length about of A ^ 1500 km. Downstream, 
1-D cuts show much stronger perturbations in the magnetic field with amplitude 
about of SB y ps O.lBo and SB Z ps O.IBq, and wave-length about of A ps 1000 km. 

The modeling gives the following range of variation in densities of the solar 
wind and pickup ions at the lunar exobase in the case with parallel magnetic 
field: n H + = 1. no, and n H + ~ (1.0 — 6.0) x 10 _4 no, n He + ~ (1.0 — 3.0) x 10 _4 no, 
n Na + ~ 1.5 x 10 _3 no- no = 3. cm -3 denotes upstream density. 

Finally, the modeling shows that the physics of mass loading is different in 
this parallel case. We find that the E-field is muted when U is parallel to B 
and that the driving E = U x B may be reduced or even disappear in this 
case. The coupling between the solar wind and pickup ions is only available 
through the excitation of the kinetic Alfven and compressional waves similar 
parallel mass loading in case of comets (see e.g. the 2.5-D modeling by Galeev 
et al. (1987)). The modeling in Galeev et al. (1987) shows the transformation 
of a ring-type velocity distribution of pickup ions with a small thermal veloc- 
ity to the shell velocity distribution. The thermal velocity of pickup ion shell 
velocity distribution becomes comparable with a local speed of the solar wind. 
The pickup ions are accelerated by the electric field perturbations excited due 
to beam-solar wind plasma instabilities in opposite the case with an oblique 
upstream magnetic field. No magnetic barrier was observed in modeling. 

The structures of plasma environment near the Moon observed in our mod- 
eling (asymmetrical Mach cone, plasma wake, pickup ion exospheres, asymmet- 
rical plasma tails and excitation of the magnetic field perturbations) need to 
be compared with the measurements by the particles and magnetic field instru- 
ments on the ARTEMIS PI spacecraft. This task can confirm the value of our 
solar wind - lunar exosphere interaction models. 
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4. Conclusions 


The hybrid modeling of lunar plasma environment with three pickup ion 
species and various orientations of the upstream magnetic field has demonstrated 
several features: 

• In the case of oblique ( 6b,u ~ 60°) and quasi-parallel (0b,u ~ 45°) up- 
stream magnetic field, the magnetic barrier has the same structure as was 
observed in the modeling in case of quasi-perpendicular upstream mag- 
netic field (Lipatov et ah, 2013b). 

• For oblique and quasi-parallel upstream magnetic field pickup ions provide 
drive formation of the asymmetrical Mach cone. The heavy Na + ions 
form an asymmetrical tail across the upstream magnetic field due to finite 
gyroradius. The light H + pickup ions form an asymmetrical tail oblique 
to the direction of upstream magnetic field. 

• In the case of parallel upstream magnetic field the mechanism of the solar 
wind mass loading with pickup ions is mostly produced by the excitation 
of the Alfven and magnetosonic waves due to beam-solar wind plasma 
instabilities. Pickup acceleration is determined by the electric field per- 
turbations from these instabilities (see e.g. (Winske et ah, 1985)). There is 
no magnetic field barrier on the sun-side lunar surface. The mass loading 
does not result in formation of Mach cone in case of parallel (0b,u ~ 0°) 
upstream magnetic field. The light and heavy pickup ions fill in the lunar 
plasma wake. 

• For upstream magnetic field in the range (0° < 0b,u < 45°), pickup ions 
generate the electromagnetic perturbations upstream and in the lunar 
plasma wake. 

• The results of this work may be important for understanding the plasma 
processes in the plasma environment near very weak comets, asteroids, 
Pluto, and Europa (Sauer et ah, 1997; Lipatov et ah, 1997, 2010a, 2013a). 
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Figure 1: Lunar plasma environment and the system of coordinates. 
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Figure 2: 0jj,b = 60°. 2-D cuts of the Na + , H + pickup ion (top and bottom) density profile. 
Here No = n H + = 3.0 cm -3 . x-y cuts (left column) are located at z = 0, y-z cuts (middle 
column) are located at x/L = 3.0, and x-z cuts (right column) are located at y = 0. 
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Figure 3: Ojj,b = 60°. 2-D cuts of the upstream (background) H + ion density (top) and B z 
magnetic field (bottom) profile, x-y cuts (left column) are located at z = 0, y-z cuts (middle 
column) are located at x/L = 3.0, and x-z cuts (right column) are located at y = 0. 
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Figure 4: 0jj,b = 45°. 2-D cuts of the Na + , H + pickup ion (top and bottom) density profile. 
Here N=n ll + = 3.0 cm -3 . x-y cuts (left column) are located at z = 0, y-z cuts (middle 
column) are located at x/L = 3.0, and x-z cuts (right column) are located at y = 0. 
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Figure 5: Ojj,b = 45°. 2-D cuts of the upstream (background) H + ion density (top) and B z 
magnetic field (bottom) profile, x-y cuts (left column) are located at z = 0, y-z cuts (middle 
column) are located at x/L = 3.0, and x-z cuts (right column) are located at y = 0. 
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Figure 6: 0jj,b ~ 0. 2-D cuts of the A/a + , pickup ion (top and bottom) density profile. 
Here No = n H + = 3.0cm -3 . x-y cuts (left column) are located at z = 0, y-z cuts (middle 
column) are located at x/L = 3.0, and x-z cuts (right column) are located at y = 0. 


18 








Y/L Y/L 



-5 0 5 10 15 

X/L 


-5 0 5 10 -5 

Z/L 


0 5 10 15 

X/L 


Figure 7: « 0. 2-D cuts of the upstream (background) H + ion density (top) and B z 

magnetic field (bottom) profile, x-y cuts (left column) are located at z = 0, y-z cuts (middle 
column) are located at x/L = 3.0, and x-z cuts (right column) are located at y = 0. 
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